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Abstract. We consider a dynamics of Higgs' field in the framework of cosmological 
models involving a scalar field •f' from Nonsymmetric Kaluza-Klein (Jordan-Thiry) 
Theory. The field & plays here a role of a quintessence field. We consider phase 
transition in cosmological models of the second and of the first order due to evolution 
of Higgs' field. We developed inflationary models including calculation of an amount 
of inflation. 



1. Introduction 

We develop in the paper cosmological consequences of the Nonsymmetric Kaluza- 
Klein (Jordan-Thiry) Theory. The theory has been extensively described in the first 
point of Ref. [1]. We refer to the book on nonsymmetric field theory for all the details 
of our theory and we present only cosmological applications of the theory. 

We discuss in the paper cosmological models involving field \P which plays a role 
of a quintessence field. We find inflationary models of the Universe and discuss a 
dynamics of the Higgs field. Higgs' field dynamics undergoes a second order phase 
transition which causes a phase transition in an evolution of the Universe. This ends 
an inflationary epoch and changes an evolution of the field ^ . Afterwards we consider 
the field ^ as a quintessence field building some cosmological models with a quintessence 
and even with a -fT-essence. A dynamics of a Higgs field in several approximations gives 
us an amount of an inflation. We consider also a fluctuation spectrum of primordial 
fluctuations caused by a Higgs field and a speed up factor of an evolution. We speculate 
on a future of the Universe based on our simple model with a special behaviour of a 
quintessence. 

The paper is divided into three sections. In the second section we discuss infla- 
tionary models and phase transitions obtained due to dynamics of Higgs' fields. In the 
third section we present an evolution of Higgs field and quintessential models of the 
Universe. We consider several approaches to get an amount of inflation. Eventually 
we calculate an amount of an inflation in a simplified model of Higgs' field evolution 
and a power spectrum of primordial fluctuations caused by a Higgs field. 
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2. Inflationary models and phase transitions 

According to new observational data [2] concerning distances of type la supernovae 
it seems that we need some kind of "dark energy" which drives the evolution of the 
Universe. This "dark energy" can be considered as a cosmological constant or more 
general as cosmological terms in field equations (in the lagrangian) . In the case of cos- 
mological model this type of "dark energy" — "vacuum energy" is a cause to accelerate 
the evolution of the Universe (a scale factor R(t)) (see Ref. [3]). The cosmological 
constant is negligible on the level of the Solar System and on the level of the Galaxy. 
Moreover it can be important if we consider even nonrelativistic movement of galaxies 
in a cluster of galaxies (see Ref. [4]). In some papers considered cosmological terms 
result in changing with time of a cosmological constant (see Ref. [5]). Some of them 
introduce additional scalar field (or fields) in order to give a field-theoretical description 
of such an evolution of a cosmological "constant" . Those scalar fields are independent 
in general of the additional scalar fields in inflationary models. 

Thus the inflation field (or fields in multicomponent inflation, which can be 
the same as some of Higgs' fields from G.U.T.-models) can be different from those 
fields. In particular considering scalar-tensor theories of gravitation results in so called 
quintessence models (see Ref. [6]). Moreover in such theories there is a natural field- 
theoretical background for an inconstant "gravitational constant" . In such a way this 
quintessence field can be used in twofold ways. First as a source of change in space and 
time of a gravitational constant. Secondly as a source of cosmological terms leading 
to the model of quintessence and a change in time of a cosmological "constant". In 
our theory we have a natural occurrence of these phenomena due to the scalar field ^ 
(or p). Let us consider the lagrangian of our theory paying a special attention to the 
part involving the scalar field i.e.: 




£kin( V \P) 



gauge 



(2.1) 




£ scal (<F) = (Mg^ + n 2 g^g 5 ^) 9 V 9 



We put c = h = 1. 

Now we rewrite the lagrangian (2.1) in the following form. 



L = R(W) - 87rG N e 
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matter 
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where in L matte r we include all the terms from Eq. (2.1) with Yang-Mills' fields, Higgs' 
fields, their interactions and coupling to the scalar field \P. 
The effective gravitational constant is defined by 

Geff = GVe-(" +2 ^ 

in such a way that the lagrangian of the Yang-Mills field in L mat ter is without any 
factor involving scalar field $ '. If the scalar field IF is constant (e.g. \P = 0) we can 
redefine all the fields in such a way that we get ordinary (standard) lagrangians for 
these fields. 

Let us consider the situation after a spontaneous symmetry breaking and simplify 
to the case of g^ = 0. We get 

L = R - 8vrG eff L matter - 8irG N L SCSll (V) + 8irG N U(V) (2.2) 

where 

£ S cai = Mgi v 9 v ■ !P 7 , M > (2.3) 



S^G N Um = -^f^lle {n - 2 ^ + e^ n+2 ^ ( ] + =e n * . (2.1) 

r \ t pl / r 



We get the following equations 

~ ^ matter seal 

- -Rg^v = 87rG eff + 8itG n T^, (2.5) 

where R^ and R are a Ricci tensor and a scalar curvature for a Riemannian geometry 
generated by g^ v = g (piy) , 

matter 

T = (p + p)u V - p^" (2.6) 
is an energy-momentum tensor for a matter considered as a radiation plus a dust. 

seal /Tf \ _ 

SvrGW = SvrGjv ( y <7^ ■ (s Q/3 ^a ■ - • «P„J + g^X cK (2.7) 

where 

X cK = 2e^*^-ll 1 V($? lt ) - e^^^R(r) - e"*-±P, (2.8) 

where mj is a scale of a mass for massive Yang-Mills' fields (after a spontaneous 
symmetry breaking for a true vacuum case). It is convenient to write 

_ e (n-2)<? e nf e (n+2)<? 

XcK = ~^—»K - —7 — P (2-9) 



(3 = /3{£,a s ,m~) 
7 = l(C,m^, a s ) 
K = 0, 1 corresponds to true and false vacuum case, i.e. 

F« t ) = 0, V(d>l t ) + 0, 

ao = 0, ai / 0. 

For a scalar field ^ we have the following equation: 



(2.9a) 



+ (n - 2) e ("- 2 )*'a x - (n + 2) e (" +2 ^/3 
- ne n *-f - (n + 2)G eff T = 



(2.10) 
(2.10a) 

(2.11) 



where T = p — 3p is a trace of an energy-momentum tensor for a matter field. For we 
are interested in cosmological models we take for a metric tensor a Robertson- Walker 
metric: 



ds 2 = dt 2 - R 2 (t) 



dr 2 



1 — kr 2 



+ r 2 d6 2 + r 2 sin 2 Od^ 2 



, fc = -l,0,l 



(2.12) 



and we suppose that I?', p,p are functions of t only. 
One gets 
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(2.13) 
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In this case we get standard equations for a cosmological model adapted to our 
theory 



3R 
~R 



1 



8^G eff -(p + 3p) - 87rGjv -(p* + 3p 



1 



1 



RR + 2R + 2k = 8TrG eff l^-(p - p) ) R 2 + 8irG N \^-{ P v - p*) ) R 2 . 

Using Eqs (2.14-15) one easily gets 

2 , , 87rG e ff 87tGat 
-H +fc = — - — pH — p^ 



(2.14) 



(2.15) 



(2.16) 



R 

where H = — is a Hubble constant 
R 

8itG nP v = 8ttG n ^-& 2 + \\ cK 

_ 2 (2.17) 
M • 1 — 

8ttG nP v = 8kG n —V 2 - -X cK . 

Let us consider a cosmological model for a "false vacuum" case, i.e. without matter 
and only with a scalar iP and a vacuum energy V(<P l crt ) / 0. In this case p = p = 
and we get 



Let us take 
Thus we get 

and 



— •• 6MR 2 d\ cl 

2M& ml* — =■ = 

R pl <M 

■2 , 1 f 1 M - 2 1, \ „o 

3 ^2 m^j 2 y 
3i? M ■ 2 1, 



& = & 1 = const. 



<N> 



(n + 2)/3xi + n7xf - (n - 2)«i = 
i? 2 + A; = ^X cl ( Xl )R 2 

-R=2 Xcl{Xl) 

x\= e . 

One gets from Eq. (2.25) 

i?(t) = fioe^ * 

where 

is Hubble constant (really constant). From Eq. (2.23) we obtain 



Xi 



-wy + y / n 5 7 2 + 4(n 2 - 4)ai/3 
2(n + 2)/3 ' 
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2.18) 
2.19) 

2.20) 

2.21) 
2.22) 

2.23) 
2.24) 

2.25) 

(*) 

2.26) 

2.27) 
2.28) 



For Qi > 0, (3 > we get 

y / n 2 7 2 + 4(n 2 - 4)qi/3 > n|-y| 



(2.29) 



„n-2 



A c i(xi) 



(n + 2) 2 /3 

n 2 7 2 - j J n 2 7 2 - js/n 2 -/ 2 + 4(n 2 - 4)ai/3 + 2a 1 (3(n + 4) 



(2.30) 



(n + 2)^ 



n 



2 7 2 - 7 V / n 2 7 2 + 4(n 2 - 4)ai/3 + 2 ai f3{n + 4) 



(2.31) 



Thus we get an exponential expansion of the Universe. Using Eq. (2.24) we get also 
k = (i.e. a flatness of a space). 

In this way we get de Sitter model of the Universe. This is of course a very special 
solution to the Eqs (2.18-20) with very special initial conditions 



R(0) 


— Ro 
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= H R 


*(0) 


= <Pi 








= 



(2.32) 



Let us disturb the solution by a small perturbation and examine its stability. Let 

R = R + 5R, \SR\,\SR\ <C R. 
One gets in a linear approximation 



2M(p + 6MH ^ + m^-j^-i^p = 



35R 



-Mip 2 



2R5R + k = ^M^ 2 (R 2 + RSR). 



One gets 



ip = ip e % Hot sin{ yjL t + 6 



Vp 4 

2 



(2.34) 
(2.35) 
(2.36) 

(2.37) 



M 

where (we take for simplicity M = 

pi 



m 



„n-2 



P 



3Mn 7 2 - 3M 7v / n 2 7 2 + 4(n 2 - A) ai (3 



2/3(n + 2) 2 M L 

+ QM ai (3{n + 4) + 4(n + Af/3 2 x\ + 4(n + 2) 2 n 2 7 x 2 - 4(n - 4) 2 ai /3 



„"-2 



7 2 /3 (3Mn + 2n 2 (n + 2) 2 ) 



2/3 2 (n + 2) 2 M 
+ 2ai/3 2 (3M(n + 4) + 2(n + 2) 

- (3M 7 /3 + 2(n + 2) 2 /3n 7 - 2n 2 7 ) y^T 2 + 4(n 2 - 4) ai /3 - 2n 2 7 2 



2/ 2 

n 



4)) 



if 



> 
(2.38) 



M > M 



2 (n 2 (n + 2) 2 7 2 /3 + 2a x (n + 2) 2 (n 2 - 4) - n 3 7 2 + n 2 7v / n 2 7 2 + 4(n 2 - 4)ai/?) 
/3 (-3 7 2 n - 6(n + 4) + 3 7v / n 2 7 2 + 4(n 2 - 4)ai/?) 

(2.39) 

Thus we get an exponential decay of the solution tp, i.e. a damped oscillation 
around = W\. It means the solution ^ = ^ is stable against small perturbations of 
initial conditions for i.e. 

W(0) = Oq + V9 sin5 (2.40) 



^(0) = 99 (-fi?osin<5 + cos5) . 



(2.41) 



Making ip and S sufficiently small we can achieve smallness of perturbations of initial 
conditions. The exponential decay of the solution to Eq. (2.34) can be satisfied also 
for some different conditions than Eq. (2.39) if we consider aperiodic case. However, 
we do not discuss it here. 
Thus we get 

< -Mip 2 < a 2 e~ 3Hot (2.42) 



where a 2 is a constant. 

From Eq. (2.35) one gets 



< SR < a 2 e- 4H °K 



(2.43) 



which means that SR ~ and SR = O(SR(0)). It means SR = 5R(0) + SR(0)t. 
Moreover from Eq. (2.36) we get that k = due to the fact that SR is a linear function 
of t. In this way we get 

R(t) * R(t) + SR(0) + SR(0)t = (fy + 6 JM±^M)l ) j {2M) 
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This means that the small perturbations for initial conditions for R result in a 
perturbation of Rq in such a way that R$ is perturbed by a quickly decaying function. 
Thus our de Sitter solution is stable under small perturbations and is an attractor 
for any small perturbed initial data. One can think that this evolution will continue 
forever. However, we should remember that we are in a "false" vacuum regime and 
this configuration of Higgs' fields is unstable. There is a stable configuration — a "true" 
vacuum case for which ax = (K = 0). For a = the considered de Sitter evolution 
cannot be continued. Thus we should take under consideration a second order phase 
transition in the configuration of Higgs' fields from metastable state to stable state, 
from "false" vacuum to "true" vacuum, from <P\ xt to ^[? rt . In this case the Higgs fields 
play a role of an order parameter. Let us suppose that the scale time of this phase 

transition is small in comparison to — and let it take place locally. We suppose locally 

a conservation of a density of an energy. Thus 

seal seal matter 

T 00 = Too + T 00 . (2.45) 

We suppose also that the scalar field & will be closed to the new equilibrium (new 
minimum for cosmological terms) and that the matter will consist of a radiation only: 

"Too = ^Ad(Ai) (2.46a) 

seal matter I q 

T 00 + T 00 = -X c0 (x ) + 8ttG n -^. (2.46b) 

For xo = is a new equilibrium point we have: 
dX C Q 



(%) = (2.47) 
= -e"* ((n + 2)/3e 2 * + n 7 ) . (2.49) 



One gets 

and supposing 
thus 




(2.51) 

/?§ (n + 2)f ~ 



(2.52) 
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From Eqs (2.46ab) one gets 



16vrG 



( x n-2 r 

n { {n + iyp l 



„_ . , ri + 2 n + 2 (2.53) 

2M 2 + I n 2 7 2 



fl™ / 1 r>\™-|-1 [ . . n + 2 _ri + 2 

(32 ( n + 2) 2 + 1 J ( n + 2) — /?— 

This gives us matching condition for a second order phase transition and simultaneously 
it is an initial condition for a new epoch of an evolution of the Universe plus a condition 
\P(t r ) = 0, R(t r ) = R(t r ), where t r is a time for a phase transition to occur. Notice we 
have simply ^ Thus we have to do with a discontinuity for a field $ '. Let us 
consider Hubble's constants for both phases of the Universe 



H 



A c i(xi) 



2 _ 

~~ 6 



H 



2 



Aco(xq) _ | 7 | 2+1 (2.54) 



Summing up we get 



1 " '" 6 3/?f(n + 2)§+ 1 

Hi + Hi 

R(t r ) = R(t r ) = R e +HoU 

#1 / ^0 
i7 2 / I 



(2.55) 



Thus we see that the second order phase transition in the configuration of Higgs' 
fields results in the first order phase for an evolution of the Universe. We get disconti- 
nuity for Hubble constants and values of scalar field before and after phase transition. 

Let us calculate deceleration parameters before and after phase transition 

RR 

5 = -fcT (2-56) 
q = -l (2.57) 

q = -1 (2.58) 



before phase transition and 



after phase transition. 

Let us come back to the Eqs (2.14), (2.16), (2.13). Let us rewrite Eq. (2.13) in 
the following form 

6MR • 2 d\ c0 



M* + + m 2 pl -^ = 0. (2.59) 
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(Remember we now have to do with a radiation for which the trace of an energy- 
momentum tensor is zero.) 

Let us suppose that & « 0. Thus Eq. (2.59) simplifies 

J_ M ^ + ^=0 (2.60) 
m^j (Mr 

and we get the first integral of motion 

-^-y 2 + A c0 = 6 = const. (2.61) 
2m 2 pl 

MV 2 = ml^S - 2A c0 ). (2.61a) 
One gets from Eqs (2.14) and (2.16) 

^ = -^•^^-2^ - Xc0 ) (2 - 62) 



Using Eq (2.61a) we get: 



3R 1 p r -= 5 , . , . 

One can derive from Eqs (2.64-65) 

^-(RR) = -\~5R 2 + l\ c0 R 2 - k (2.66) 
at 6 6 

A c0 = -e^(/3 2 *+ 7 ). (2.67) 

Let us take k = and (5 = and let us change independent and dependent variables 
in (2.66) using (2.67) and (2.61a). One gets 

2y 2 (y 2 - 1)0 + y ((n + 4)y 2 - (n + l)y - 2) = ^(y 2 - l)f(y) (2.68) 



where 



(2 . 69 ) 
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and 

It is easy to see that y 2 > 1 for 



f = R . 



(2.70) 



(2.71) 



Moreover according to our assumptions & and this can be achieved only if < 
y — 1 < e, where e is sufficiently small. 

Thus for further investigations we take z = y — 1, y = z + 1. One gets 



d 2 f df 
2z(z + 2)(z + 1 f^ + ( z + !) ((" + 4 )^ 2 + ( n + 7 )^ " (« + 3 )) ^ 

4M 

-_s( z + 2)/(z) = 0. 



(2.72) 



Let us consider Eqs (2.64-65) supposing k = and (5 = 0. 

After eliminating f? (via differentiation of Eq. (2.65) with respect to time t) and 
changing independent and dependent variables one gets: 



d_ 

dy 



8vrG 



N 



d_ 

dy 



(f 2 )y n (y 2 - 1) 



where 



Pr = Pr 



It is convenient for further investigations to consider a parameter 

AM 



r = 



In such a way we get 

2y 2 (y 2 - 1)0 + y ((n + 4)y 2 - (n + l)y - 2) ^ - r{y 2 - l)f(y) 



and 



2z(z + 2)(z + l) 2 ^l + (z+ 1) ((n + 4)z 2 + (n + 7)z - (n + 3)) f- 

az z ' dz 

-rz{z + 2)f(z) = 0. 
One can transform Eq. (2.73) into 



(2.73) 
(2.74) 

(2.75) 

(2.76) 

(2.77) 



d_ 

dy 



8nG- Pr 1 3 ""'~ 2 



,n+2 2 



+ - 1) / 



y™- 1 [(n + 2)y 2 -n]/ 2 . (2.78) 
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Changing the independent variable from y to z = y — 1 one gets: 



d_ 

dz 



^ G ^ z+ Hr i)n+ 2 + 2 {Z+1)nZiZ + 2) ) f 

= (z + l) n - 1 [(n + 2)(z + l) 2 -n] f. 



(2.79) 



Let us come back to Eq. (2.71) in order to find time-dependence of \P. One gets 



±- 



t-u 



or 



/ 



dx t — £i 

VM ' 



x^ f3x n+2 — \^f\x r 



(2.80) 
(2.80a) 



If n = 21, where I is a natural number, one gets for (3 > 

2p+l 



M )""(EE C)( " 2)fc " P 

1/ I p=0 



2p 



+ 



1 J j^jn^3 )| _ ± (*-*i) 



(2.81) 



2^3 V\/^-|7l + V3(^ + l7l) 



VM 



In this case it can be expressed in terms of elementary functions. For a general 
case of n, (3 and 7 one gets 



2*-f 1 F 2 (i,-f,(l- f);^) _^( f _ fl) 



V7™ 



VM 



(2.82) 



where 1^2(0, &, c; z) is a hypergeometric function. For a small z (around zero) one finds 
the following solution to Eq. (2.77) 



r 



(n-1) 



(2.83) 



Ci,C2 = const. Thus 



R(y) 



c,fa-i)^(i- ( " + J )( ^t 23) (y-D 



n + 1 



+ C 2 1 + 



(n-1) 



(2.84) 
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where y > 1 (but only a little). 

Let us make some simplifications of the formulae taking under consideration that 
y — 1 is very small. 

One gets 



f(y) = C (1 + 

R(y) =Rt(l + 
After some calculations one gets 
8wG N p r = y 2 ^ 



(n-1) 



(n-1) 



(y-1) 2 , C>0 



(2.85) 
(2.86) 



2(n + 5) ^ 4 + 10r(ra + 2) ^ 3 



(n + 4)(n-l) (n-l)(n + 4)' 
(8r- in- l)(n + 4)) 9 2r(n + 3) (r-n + 1) 



2(n- l)(n + 4) 



(n 2 - 1) 



(n-1) 



(2.87) 



(taking into account that (y — 1) is very small). Let us make some simplifications in 
the formula (2.71) for y closed to 1. We find: 



i i rc+ 4 / — 

_^ = = ± JlLL^ (t _ tl) 



and finally 



■ 71 + 4 



y 



■(t-ti) + l 



(2- 



(2.89) 



such that y = 1 for t = t\. 

Let us pass to the minimum value for p r with respect to y closed to 1. Thus we 
are looking for a minimum for a function 



1 



8vrG 



-y 



2(n+l) 



2(n + 5) 4 10r(n + 2) s 

ry + t — TT7 — 



L (n + 4)(n-l) w (n-l)(n + 4)' 
(8f- (n- l)(n + 4)) 2 2r(n + 3) (r-n + 1) 



2(n- l)(n + 4) 



(n 2 - 1) y (n-1) 



8ttGn 



y 



2(n+l) 



W 4 (y) 



(2.90) 

for such a y that is closed to 1. 

Let us write y = 1 + n, where is very small and develop W4(l + n) up to the 
second order in 77. One finds 

W 4 (l + n) ~ V 2 (rj) (2.91) 



and 



V 2 (n) = an 2 + 6n + c 



(2.92) 
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where 

_ (-n 2 + 3n(20r-9)+4(32r-31)) 

a ~ 2(n + 4)(n-l) [ j 

_ (-n 3 - 4n 2 (7r - 3) + 5n(18r - 11) + 22(3r - 2)) 

b = (n + 4)(n2-l) (2 ' 94) 

_ (-3n 3 + 2n 2 (9r - 8) + n(50r - 29) + 4(7r - 4)) 

C= 2(n + 4)(n 2 -l) (2>95) 

Let us find minimum of ^(77) for 77 > such that V^r?) > 0. For a < 0, b < 0, 
c > we get 

F 2 (??mi„) = (2.96) 

and 

= (2.97, 
2\a\ 

In this case 

(1 + r/ min ) 2(n+1) ^4(l + r/ min ) > (2.98) 

and will be closed to the real minimum of the function (2.90). Simultaneously 

??min < 1- 

For 

Pr = &oT* (2.99) 



(where a is a Stefan-Boltzmann constant) we get T m ; n and we call it a Td (a decoupling 
temperature-decoupling of matter and radiation). Thus yd = 1 + r? m in- Ud is reached 
at a time td (a decoupling time) 

„ ri 

" = *> + *nX- ( 2 ' 10 °) 

2 ry K 2 J 

From that time (td) the evolution of the Universe will be driven by a matter (with 
good approximation a dust matter) and the scalar field \P. The radius of the Universe 
at t = td is equal to 

R(t d ) = i?l (l + ^^min) 2 • (2-101) 

The full field equations (generalized Friedmann equations) are as follows: 

(r) ^^ + 3 2^ + 2H < 2 - 102 > 



3R _ SttGW \[~M_ x j/i 
i? " 3 Pm 2 m 2 



Pm - - ( 2—^ - A c0 I (2.103) 



2M ~ 6M R- d\ c0 . 

—V ~ — ■ -5* ~ + n + 2 8vrp m G eff = 0. 2.104 
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Before we pass to those equations we answer the question what is a mass of the 
scalar field & during the de Sitter phase and the radiation era. One gets: 



m? = ^^(ff!) (2.105) 

1 2M d& 2 y ' y ' 

2 m pi d 2 X c0 . . 

m ° = mn^ {%) - (2 - 106) 

The second question we answer is an evolution of the Universe during a period from 
t r to t\, i.e. when ^ = & , i.e. y = to y = 1. 

This will be simply the de Sitter evolution 

R(t) = Re Hlt (2.107) 

where 

Hl = Jh^A . (2.108) 

This will be unstable evolution and will end at t = ti, i.e. for \P ~ and value of the 
field ^ will be changed to the value corresponding to y = 1, i.e. 

nti) = \^^- v (2-109) 

Thus we have to do with two de Sitter phases of the evolution with two different Hubble 
constants Hi and H . In the third phase (a radiation era) the radius of the Universe 
is given by the formula 

1+ _ 171 (t-h) 2 ) (2.110) 

3M(n- l)(3 n J 

and such an evolution ends for t = tj, 



r 



R(t d ) = R 1 exp(H t r + H 1 (t 1 -t r ))-^l + -^—^r 1 2 min j . (2.111) 

In both de Sitter phases the equation of state for the matter described by the scalar 
field & is the same: 

pq) = -p#. (2.112) 

In the radiation era we get 

p* = ip*. (2.113) 
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Now we go to the matter dominated Universe and we change the notation intro- 
ducing a scalar field Q such that 

Q = v / M— . (2.114) 



In this case we have 



m p i 



PQ = \Q 2 + U{Q) (2.115a) 
PQ = \Q 2 ~ U{Q) (2.115b) 



where 



( n _|_ 2) n 
a = — m p i, 6 = — = m p i. (2.117) 
VM VM 

Let us write Eqs (2.102-104) in terms of Q: 

Q + 3i7Q + U'(Q) - 8iT(n + 2)G eff/ 9 m = (2.118) 



where 



(ra + 2) „ 



G eS = G N e P (2.119) 

H 2 = 1 (^87rG eff/ 9 m + ^Q 2 + [/) (2.120) 

^ = 1 (8vrG effPm + 2(Q 2 - 17)) . (2.121) 
Let us define an equation of state 

and a variable 

_1 + Wq_IQ 2 

XQ - x^W Q - IT' (2 ' 123) 

which is a ratio of a kinetic energy to a potential energy density for Q. It is interesting 
to combine (2.118) and (2.120) using (2.122) and (2.123). 
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Let us define 



±o 2 + U 

°Q = 4^ (2-124) 



3H 2 



87rG eS p m 

= 3H 2 — • (2.125) 
One gets 

Q m + Qq = 1 (2.126) 

or 

8nG eS p m = 3(1 - Q Q )H 2 . (2.127) 

We call p m an energy density of a background and for such a matter the equation of 
state is 

W m = ^ = 0. (2.128) 

Pm 

It is natural to define 

p m = 87rG eS p m (2.129) 

as an effective energy density of a background with the same equation of a state 
as (2.128). 

Under an assumption of slow roll for Q, i.e. 

X -Q 2 « U(Q) (2.130) 

one gets 

W Q ~ -1, (2.131) 

i.e. an equation of state of a cosmological constant evolving in time. In this case 
xq ~ 0. This can give in principle an account for an acceleration of an evolution of 
the Universe if we suppose those conditions for our contemporary epoch. 



3. Evolution of Higgs Field and Quintessential-Cosmological Models 

Let us come back to the inflationary era in our model. For our Higgs' field is 
multicomponent, i.e. it is a multiplet of Higgs' fields, we have to do with so called 
multicomponent inflation (see Ref. [7]). In this case we can define slow-roll parameter 
for our model in equations for Higgs' fields. 

One gets from the first point of Ref. [1] (see Eq. (5.7.7), p. 385) 

i(f).-«"( I f). = -T^{s*} («) 

^ n ) av 
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where Lq~ = L d d (a time component of L d ~) is denned by 



V, — F and I \ are defined in [1] (Eq. (5.6.9), (5.6.10), (5.424)). 

h y n ■> an 



/ (n-2)* 2e- 2 * • 1 \ 

3H 2 = 8ttGn {—pr- V '(*) ~ -p-4d=(*) + 2 A C o(^o)J 



2e~ 2q, ° 

H = 8itG n ( — £ ki „(<£) 



where H is a Hubble constant and 



ob dt' 

<?o is a constant value for a field 1^. The slow-roll parameters are defined by 



P _ ; - av 

H 2 

5 = 



lab \ 


1 n bfi 


d | 


&)*«)}.. 


H 2 i 


[lab{ 


[ 9 i8 %i (*a}j 



with the assumptions of the slow roll 

e = 0( V ), 5 = 0( V ) 



for some small parameter rj. 

Under slow-roll conditions we have 



P n& f at/' \ 



n / av 



8tt ( p( n - 2 )^o ~] 1 



and with a standard extra assumption 

8 



jj = 0{r, 2 ) 
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We remind that the scalar field \P is the more important agent to drive the inflation. 
However, the full amount of time the inflation takes place depends on an evolution of 
Higgs' fields (under slow-roll approximation). The evolution must be slow and starting 
for = d>l rt . It ends at <P = <P° rt . 

Thus we have according to Eq. (2.26) 

^ = ^ n p/i 6nd \ = H (t en d — initial) 
^HlinitialJ 



where 

<P(t end ) = <2> c ° rt (3.12) 

^initial) = ^ert (3-13) 
^end t r 

(we omit indices for <P). 

Thus we should solve Eq. (3.9) for initial condition (3.12) and with an assumption 
H = H . We have of course a short period of an inflation with H = Hi from i en d = t r 
to t = t\, i.e. 

JV tot = N + N 1 = H (t r - initial) + " *r)- (3.14) 

Let us remind to the reader that N is called an amount of inflation. Let us consider 
inflationary fluctuations of Higgs' fields in our theory. We are ignoring gravitational 
backreaction for the cosmological evolution is driven by the scalar field & '. We write 
the Higgs field as a sum 

$% l (r,t) = $° h (t)+5<I><? h (r,t) (3.15) 

where ^^(t) is a solution of Eq. (3.9) with boundary condition (3.12-13) and 5^^(r, t) 
is a small fluctuation which will be written in terms of Fourier modes 

K 

The full field equation for Higgs' field in the de Sitter background linearized for S$ a gAt) 
can be written 



4i(mQ +3H ±(M±r) +1-K 2 (m'C] 

dt 2 \ Kb )av >dt\ Kb )av R 2 V KJ 

' 5 2 V (#?(i) 



2r 



9if. 



(3.17) 



where 



l dc M d Ri + l cd9 arng™ c Mi d = 2l cd g drh g™ '6$%.. (3.18) 
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Let us change dependent and independent variables. 

dr=^-, R(t) = Roe Hot (3.19) 
K(t) 

JKr) = --L (3.20) 

— oo < r < (a conformal time) (3-21) 

and 

a 

5$% . = -JgL . (3.22) 

Km Jl v ' 

Simultaneously we neglect the term with the second derivative of the potential. Even- 
tually we get 

d 2 (~ A \ 2 d 



dr 2 



where 



l dc M K ~ b + l cd9 amg mc M a g . = 2l cd g^g mc X i £ . (3.24) 
Let us notice that \Kt\ is the ratio of the proper wavenumber \K\R~ 1 to the 

Hubble radius — . At early times \Kt\ S> 1, the wavelength is small in comparison to 
Ho 

Hubble radius and the mode oscillates as in Minkowski space. However, if r goes to 
zero, \Kt\ goes to zero too. The wavelength of the mode is stretched far beyond the 
Hubble radius. It means the mode freezes. Thus the analyzes of the modes are similar 
to those in classical symmetric inflationary perturbation theory neglecting the fact that 
we have to do with multicomponent inflation and that the structure of representation 
space of Higgs' fields should be taken into account. 

However, the form of Eq. (3.23) is exactly the same as in one component symmetric 
theory if we take a field to be considered. The relation (3.24) between Mi and 



y - is linear and under some simple conditions unambigous 

Ka 



d 

a 

It is possible to find an exact solution to Eq. (3.23) and Eq. (3.24). One gets 

t\K\ cos(r|i?|) — sin(r|i? : |) \ 



and 



yl = C a - 
A Km lKm 



+ C 2Km 



Mi- = N d ~~ 

Kb 1Kb 



~ 2Kb 





T 




't\K\ 


sin(r|i?|) 


+ cos (ri?) 




r 




t\K\ 


cos (ri?) 


— sin (ri?) 




r 




t\K\ 


sin(r \K |) 


+ cos (ri?) 



(*) 



(**) 
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where 



l dc Nf Rl + 1*9**9™ N^. = 2l cd g arh g™ £ C? R „, i = 1,2. 



In this way one obtains 



(***) 



Km 



Rq 



K ™\ r h 



(****) 



and using (*) one easily gets 
( \K\ 



5$ 



C" 



Km lKm \ 

\K\ 



e~ Hot cos 



\K\ 

RqHq 



-Hot 



+ C 2Km 



\K\ 



RqH 



■ sin 



RqH 



-H t 



+ COS 



+ sin j 

\K\ 
RqHq 



\K\ 



-Hot 



RqHq 
-Hot] 



(V*) 



where 



(j a ^ (J a _ 

iKm " iKm 



(VI*) 



However, it is not necessary to use a full exact solution (V*) to proceed an analysis 
which we gave before. 

Thus the primordial value of R K {t) is equal to 



R k ~ ~ 



H ° -60%. 



d_(<pd \ Kr 

dt y^m) 



\t=t* 



(3.25) 



where R K (t) is defined in cosmological perturbation theory as 

R(V(t)=4^R l} (t) 
(see Ref. [8]) and t* is such that \K\ = R(t*)H , i.e. 



(3.26) 



t* = — In 
Mo 



\K\ 



RqH 



(3.27) 



The primordial value is of course time- independent. We suppose that fluctuations of 
Higgs' fields (the initial values of them) are independent and gaussian. 

Thus we can repeat some classical results concerning a power spectrum of primor- 
dial perturbations. One gets 



\ U Krh\ 



Ho 
2vr 



(3.28) 
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Thus 




t=t* 



for t* given by Eq. (3.27). 

Taking a specific situation with concrete groups H, G, Gq and constants £, £, r we 
can in principle calculate exactly Pr(K) and a power index n s . Thus using some 
specific software packages (i.e. CMBFAST code) we can obtain theoretical curves of 
CMB (Cosmic Microwave Background) anisotropy including polarization effects. 

Let us give the following remark. The condition for slow-roll evolution for tf> can be 
expressed in terms of the potential V. If those are impossible to be satisfied for some 
models (depending on G, G ,H, G' and parameters £, £), we can employ a scenario with 
a tunnel effect from "false" to "true" vacuum and bubbles coalescence. In the last case 
the time of an inflation (in the first de Sitter phase) depends on the characteristic time 
of the coalescence. 

Let us consider a more general model of the Universe filled with ordinary (dust) 
matter, a radiation and a quintessence. 
From Bianchi identity we have: 

d 3R . . 

^Ptot + (Ptot +Ptot)-^r = (3.30) 

where 

Ptot = PQ+Pm+Pr (3.31) 

Ptot=PQ+Pr- (3-32) 

Let us suppose that the evolution of radiation, ordinary matter (barionic + cold dark 
matter) and a quintessence are independent. In this way we get independent conser- 
vation laws 

d (R 3 p m ) = (3.33) 



dt 
d_ 

dt 



d (R%) = (3.34) 



and 



One gets 



PQ + (PQ+PQ)^ = 0. (3.35) 



^4 

Pm = ' A = const ' ( 3 ' 36 ) 

p r = — , B = const. (3.37) 
R 
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From Eq. (3.35) we get 



p Q + (l + W Q ) PQ -^ = 0. (3.38) 



Substituting Eqs (3.36-37) into Eqs (2.14) and (2.16), remembering Eqs (2.114) and 
(2.116-117) one gets 



R.\ . , I (A B \ I 
R 



+ k = 


1 




3m 2 pl 


3R 


1 


~R 





and from Eq. (2.118) 



Q + 3HQ- = 0. (3.41) 



However, now Q and U (Q) are redefined in such a way that — =- factor appears be- 

m 2 p] 

fore pq. This is a simple rescaling. We have of course 

p m = e- {n ^Pm = e~ aQ Pm (3.42) 

p r = e -( n+2 ^p r = e~ aQ p r . (3.43) 

Thus we get 

Pm = ^e a Q (3.44) 

Pr = J^°«. (3.45) 

Let us consider Eqs (3.39-41) and Eq. (3.35). Combining these equations we get 
the following formula 

(n + 2) A 



m pl R 3 



(1 + W Q )p Q = 0. (3.46) 



The one way to satisfy this equation is to put Wq = — 1. It is easy to see that if A = 
then Eq. (3.46) is satisfied trivially. The condition A = does not imply B = and 

VM 

we can still have B / 0. For such a solution Q = Q = const, and Qo = \Pq found 

m p i 

earlier 



X ° \j{n + 2)f3 a s \m pl ) \j ( n + 2)R{r)' (3 ' 47) 
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Thus 

n + 2 



n+2 



n + 2 



e - aQo = e -(n + m = a n + 2 (^y ^ ^ 2 imp . (348) 



or 



e" aQo = < +2 ( f- ) ( — ) 1^1 • (3-49) 



i)(^r(¥) 

Let us put Wq = -1 into Eqs (3.39-40). One gets 

1 (A B 



n + 2 
2 



:!/.' 1 / A 2B 

+ ^M+PQ)- ( 3 - 51 ) 



R " \2R 3 3R 4 

From Eq. (3.50) one gets 

±RdR 



= dt (3.52) 

\j 3^"^ ~ + 3^ R + 3^% 

if A = 0, 



= cft. (3.53) 



o 2~~ — fc-R "I - o T" 

3m , 3m , 

P l P l 



Taking x = R 2 one gets 

/ , pQ ^ = ±2(/ - ^o) (3.54) 

and finally 

/ , * = ^»( t - t „) (3.55) 



where 




x=\l — y. (3.56) 
PQ 



For a flat case k = we find: 

_ 2\/3^/p^mpi 
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(t-to). (3.57) 



The last formula can be easily integrated and one gets 



f>' — i 4 / — A 



(3.58) 



where we take a sign + in the integral on the left hand side of Eq. (3.57). 
Let us calculate a Hubble parameter (constant) for our model. One gets 



H = | = V3m pl (^p) ctgh 



B 



(t-t ) . 



(3.59) 



For a deceleration parameter one obtains 

RR (2sh 2 (u) - ch(u)) 



q = - 



R 2 



ch(it) sh 2 (it) 



where 



2VV7^m pl 
u= w — (t-t ). 



(3.60) 



(3.61) 



Thus we get a spatially flat model of a Universe dominated by a radiation which is 
expanding and an expansion is accelerating. 

Let us take k = and B = in Eq. (3.52). One gets 



±RdR 



I PQ 

V K, 



eft 



i? 4 + 



3 m 



i? 



pl 



PQ V ^4 



One finally gets 



Let us notice that 



/ 



x dx 



V^x 3 + 1) 3m p] 



(3.62) 
(3.63) 

(3.64) 



R 
R 



A 



PQ 



RR _ (2p Q R 4 - A)R 
~W 2(A + PQ R3) 

Thus the model of the Universe is expanding and accelerating. 



(3.65) 
(3.66) 
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Let us consider the integral on the left hand side of Eq. (3.64). One gets after 
some tedious algebra: 



^| ^92^4\/3 + 3 (l41\/3 + 272) + ^9>/3 - 15) 16 + 9Vsj F 



(3.67) 



where 



W 



2 2 Ua:-V3-1 



+ 2\/6-\/3 



'(2g + V3-1)' 
(2x - \/3- 1) 



+ 7-4\/3 



(2a + \/3-l) 
(2a; - y/3 - 1) 



1 + 



V3 



II = II arccos 



(2x + V3-l) \ 
(2s- V3-1)/ 

' (2x + (y/3- 1)) 

i (2a;-(V3+l)) < 

V (2s + (V3-l)) ' 
, I (2a:-(V3 + l)) 



(3.68) 



2(2 + V3) ,(V3-l), 



(5 - 2^3) 



13 



F = F \ arccos 



2(2 + V3) , 



(5 - 2>/3) 
13 - 



(3.69) 
(3.70) 



where FI is an elliptic integral of the third kind and F is an elliptic integral of the first 
kind 



u 



dp 



y/l - k 2 si 



sin 2 ip 



(5 - 2^3) 



13 



n(u 



,n,k) = / — 
J ( 1 



dp 



(1 — nsin 2 <p)-\/l — k 2 sin 2 99 



n = 





V3-1 



(3.71) 

(3.72) 

(3.73) 
(3.74) 
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and with the same modulus as F(u, k) 

k 



o 5- 2^/3 



13 



However, the most important condition for an existence of a physical solution is coming 
from the first part of a sum in the right hand side of Eq. (3.67). 
We should have 

W > 0. (3.75) 
In order to analyze condition (3.75) let us write W in the following way 



_ 20 _ 9^1 + (16^1)^ + 2 ^T7^ {y + 7 _ 4^3) (y-l + 4) 
W = " (3.76) 

y- 1 

where 

Obviously y > 0. 

We have two possibilities: 

I y>\ (3.78) 

and 



-20 - h£ + ^~l^ y + 2V 6 - v^y (y + 7 - 4^3) (y - 1 + ^) >0 (3.79) 

II y < 1 (3.80) 
and 

-20 - 2yJ + (16-7^ y + 2 ^ 6 _ + 7 _ 4v /3) ^ _ i + ^ < o, (3.8i) 
The first possibility can easily be solved: 

I y > 4.0612 (3.82) 
or for the second possibility 

2 \/3 

II < y < 1. (3.83) 
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We have finally 



For 



0.7916 < x < 



v / 3+ 1 



3^-5 1 

^ < x < -. 

2 2 



x > 



we get 



x < 3.072. 



Let us notice the following: the function W has a finite limit at 
we can remove a singularity at x = v/ ^ +1 . In this way we get 

0.7916 < x < 3.072 



or 



3^-5 1 
—^ <X< 2- 
This simply means that the solution exists only for 



or 



A I A 

0.7916 4 3 /— < R < 3.072 ? 

PQ V pq 



2 V pQ 2 \l PQ 

Let us calculate a density of a quintessence energy: 

PQ = ml l U{Q ) = - l -e^{^-\ 1 \) 



n 



n + 2 



m , 



2(n+l) 



a 



P 



m p J \R(r) 



An energy density of a matter is equal to 

A 



and 



Hi ~ e aQo 



a 



n+2 



11 



n + 2 



2 
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Let us consider the behaviour of the "effective" gravitational constant in our model 



G eff = G N e~^ +2 ^ = G N e- a< 3° = G N (±) 



n+2 



■ / ~ ~ N. n + 2 



(3.95) 



2 



Thus G e ff is a constant. For we are living in that model we should rescale the constant 
and consider G e s(Qo) (Eq. (3.95)) a Newton constant. Thus 



n + 2 n + 2 ----- - H+2 



where Go is a different constant responsible for a strength of gravitational interactions 
in earlier epochs of an evolution of the Universe. It means we should write 



n + 2 n + 2 - -- - 2±2 ' 



°-= (°-' 2< " +2> (^) " (^) " (Aj )- ( " +2) * ^ 

The obtained solution (i.e. (3.67)) is for Wq = —1. However, this is not an attractor 
of the dynamical equations. This is similar to the tracker solution of Steinhardt (see 
Ref. [9]). Moreover we cannot apply his method because our equation for a scalar 
field \P (or Q) is different. It contains an additional term with a trace of the energy- 
momentum tensor (matter + radiation). Only with a radiation filled Universe our 
equation is the same (if we identify p r with Steinhardt density of a matter). In this 
case we could apply Steinhardt tracker solution method and apply his criterion for a 
potential U(Q). In our case U(Q) is given by Eq. (2.116) 

U'(Q) = — ^ Un + 2)/3e aQ - nh\e bQ ) (3.98) 
2VM 

U"(Q) = — ((n + 2) 2 f3e a ® - n 2 \ 7 \e bQ ) . (3.99) 
2v M 

The Steinhardt criterion consists in finding 



U"U _ {(n + 2) 2 f3e a ® -n 2 | 7 |e^) ((3e a ® - | 7 |e 6 «) 
" (U') 2 ~ ((n + 2)(3e a Q -n\j\e b Qf ' 



To get a tracker solution we need 

r>i. (3.101; 
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In our case 



r = 1 - 



4/3 | 7 | e (a+fe)Q 



((n + 2)/3e a <2 - n\^\e b Qy 



r ~ 1 - 



4| 7 | 



(n + 2)/3 



,(6-o)Q 



4| 7 | 



/3(n + 2) 



2m P i ) 



4| 7 | 



-2<? 



/?(n + 2) 



(3.102) 



(3.103) 



Thus asymptotically we are closed to r = 1. Moreover the Steinhardt criterion is 
only a sufficient condition. Our radiation filled Universe seems to be unstable due to 
appearing of a matter (a dust matter). Moreover a quintessential period of a Universe 
model has a severe restrictions of a value of a radius. Thus the solution with Wq = — 1 
and a matter independently evolving cannot evolve forever. 

In order to answer a question what is a further evolution of the Universe we come 
back to Einstein equations with a scalar field ^ and matter sources. Supposing as usual 
a Robertson- Walker metric and a spatial flatness of the metric we write once again a 
Bianchi identity 



' seal 



+ 



1 



i t matter 

-(n+2)<? rp fxv 







rn 



pi 



m 



pi 



(3.104) 



where 



matter 

rp t^is 



is an energy-momentum tensor for a dust matter. One gets 



(3.105) 



rn 



pi 



M 



-V 2 — = 0. 



rat 



We make the following ansatz concerning an evolution of a matter density 

Po = const. 



R 



(3.106) 



(3.107) 
(3.108) 



In that moment the scalar field and the matter p m interact nontrivially. 

Using Eqs (3.107-108), (3.106) and an equation for a scalar field & one obtains 



Po_ 

2 



m 



(n + 2) po 3R 
o — Po H — — = 0. 



pi 



m 



pi 



m li R 



(3.109) 
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Using Eq. (3.109) and Eq. (2.102) one gets 



* ± ^ (!i±I) (3.110, 
dt V3M V V m pi J m P i 2 

I , ^ = ±-^(* - to) (3.111) 



n + 1 \ 2 3p 



where 



5 = 2 — --^ (3.112) 

z = e*. (3.113) 
Using Eq. (3.111) and Eq. (2.102) one gets 

±lnR=^±V (3.114) 
dt m p \ 

or 

R = R e { ^ ){t ~ to) . (3.114a) 

Let us consider Eq. (3.111) and let us suppose that x < 1 (and practically small). In 
this case we can simplify 

r dx ^ r dx 

J x^S + 3/3x n + 2 + 3 7 x" J x^6 + 3-/x n 

and finally we get 



»-s i, (sffl)"^VS^ ( '-* ,) - (3 ' 116) 

Thus our prediction that x < 1 and is small has been justified. Using Eq. (3.116) one 
easily gets 

*=# + w exp t-^' (1 - ( » ) ) (3 ' 117) 

"-'^-^^(-{w^-'A- (3 - iis) 



It is easy to see that 
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1. (3.119) 



Thus in this model we have to do asymptotically (practically very quickly) with a stiff 
matter equation of state: 

p = p. (3.120) 

Simultaneously an energy density for a scalar field & is dominated by a kinetic 
energy which is practically constant and 



2 V " 3n 2 



(3.121) 



In this sense a dark matter generated by ^ in the model is in some sense if-essence 
(not a quintessence). 

Let us calculate an energy density for a matter 

Q + 2) / X 

Pm=Po(Jn) " exp [ -^±^-J^=m pl ( t -t ) ) . (3.122) 



3|7l 



n 



3M 



The effective gravitational constant reads 



G e ff — Gn 



3 



exp (^) vm mpl{t - to) 



(3.123) 



It is easy to write R as a function of 



R = Rq exp 



n+l /= 



M 



pi 



2^(5 + 3(3x n+2 + 3-fX % 



(3.124) 



where 



x 



3|7l 



2n 



(3.125) 



Eq. (3.124) can be easily reduced to the Eq. (3.114) using Eq. (3.111). 

Thus an energy density of matter goes to zero in an exponential way. The effective 
gravitational constant is growing exponentially. What does it mean for the future of 
the Universe? First of all the Universe will be very diluted after some time of such an 
evolution. Secondly, a relative strength of gravitational interactions will be stronger. 
This means that if we take substrat particles as cluster of galaxies then in a quite short 
time any cluster will be a lonely island in the Universe without any communication 
with other clusters. All the clusters will be beyond a horizon. On the level of a single 
cluster the strength of gravitational interactions will be stronger and eventually they 
collapse to form a black hole. In an intermediate time a gravitational dynamics on the 
level of galaxies' clusters will be governed by a Newtonian dynamics (nonrelativistic) 
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with a gravitational potential corrected by a cosmological constant. The cosmological 
constant induces a positive pressure (a stiff matter) and will be important up to a 
moment of sufficiently large gravitational constant. After a sufficiently long time only 
a classical Newtonian gravity will drive a dynamics of galaxies' clusters. Let us roughly 
estimate this effect. In order to do this let us suppose that 



= Pi> 

In an energy momentum we get 



3n 2 



const. 



(3.126) 
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However, we have in the place of Newtonian constant Gn a G e s which is a function 
of cosmological time and in the equation for gravitational field for g^ v we should put 



1 seal 
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3n 2 



{2vPu v - g»") 



(3.128) 



where 



m p i 



x exp 



m p im; 



n+2 ri + 2 

^2 (n + l) 2 -3p ' 
3|p| " 



(3.129) 



Thus 



where 



] seal 
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T = Ae- K( *-* n) (2uV - 



pi 



(3.130) 



H 

3ra 2 



n + 2 



2(5 \ 
3lF J 

-(n+2) 



2(n + l) 2 - 3p c 
3|p| 



(3.131) 
(3.132) 



Consider the following case (the justification will be given below): 



(m eff)2 



(3.133) 
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Now let us solve the following problem. Let A be a value Ae~ K ^~ to ^ in an established 
cosmological time t = t and the same for G = G e s(t). Let us consider a gravitational 
field of a point mass Mq static and spherically symmetric. Consider a metric 

ds 2 = B{r)dt 2 - - r 2 (d6 2 + sin 2 9 d9 2 ) (3.134) 

B(r) ' 

in spherical coordinates and Einstein equations 

Rn» ~ l^g^R = ~Ag^ (3.135) 

for the metric (3.134) where 

A = Ae- K{i ~ to \ (3.136) 
Via a standard procedure one obtains 

B(rH1 _^_*!. (3 , 37) 

Thus an effective nonrelativistic Newton-like gravitational potential has a shape 

„ (r) + (3 . 138 ) 

r 6 

Remembering that G is growing exponentially in time and A is decaying in time ex- 
ponentially we see that the second term is important only for a short time. For a 
contemporary time 

A ~ c 2 • l(T 52 m 2 (3.139) 

according to the Perlmutter data and it could be important on the level of a size of 
galaxies' cluster. However, in a short time it will be negligible even on that level. 
Let us calculate a Schwarzschild radius for a mass M with an effective gravitational 
constant G e ff. One gets 



r 

9 



= J^toU = ^2M G eS (t) = r g (^p) 
3|p| 



2(n + l) 2 - 3p 



n + 2 r , , v -. (3.140) 

/l fn + 2\ [25 \ 

exp 



2 V 2 / V 3M 



Thus r| ff (calculated for an effective gravitational constant G e ff(t)) is growing 
exponentially. It means that after some time it will be of an order of size of galaxies' 
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cluster. Then galaxies' clusters collapse to form black holes. The time of a collapse is 
easy to calculate. 



r g — '"cluster 



Collapse = to + (n + 2) y/0M 

'n + 2 



+ 



In 



ln|^U(n + 2)ln(^p) 

-l 



2(n+ l) 2 -3p 



V5(n + 2)m pl (3.142) 



(3.141) 



where r g is a gravitational radius of a mass of a galaxies' cluster with G = Gjy. 

Let us calculate a Hubble constant and a deceleration parameter for our model. 
One gets 



R 

H = — = (n + l)m p i 



RR 



-1. 



(3.143) 
(3.144) 



Summing up we get the following scenario of an evolution. The Universe starts in a 
"false" vacuum state and its stable evolution is driven by a "cosmological constant" 
formed from R(r), P, V(^crt) and $i being a minimal solution to an effective self- 
interacting potential for <P\ This evolution is stable against small perturbation for \P 
and R. The evolution is exponential and a Hubble constant is calculable in terms of 
our theory (the nonsymmetric Jordan-Thiry theory) . This evolution ends at a moment 
the "false" vacuum state changes into a "true" vacuum state. In that moment an en- 
ergy of a vacuum is released into radiation (and a matter). The Universe is reheating 
and a big-bang scenario begins at a very hot stage. Moreover the evolution is still 
governed by a scalar field & which attains a new minimum of an effective potential. 
The effective potential is different and a new value of a is different too. The 
evolution is exponential and a new Hubble constant Hi is also calculable in terms of 
the underlying theory. In this background we have an evolution of a multiplet of scalar 
fields from <£* rt to #° rt . This goes to a spectrum of fluctuations calculable in the 
theory. After that time the field \P is slowly changing iP ~ and the temperature of 
the Universe is going down. The evolution of a radiation is in the second inflation- 
ary stage adiabatic and afterwards during next phase (p ~ 0) nonadiabatic. The last 
means there is an energy exchange between radiation and a scalar field The total 
amount of an inflation iVtot is calculable in the theory. The evolution of a factor R(t) 
is governed by a simple elementary function of a small rise. After this the radiation 
condenses into a matter (a dust with zero pressure) and a matter (a dust), a radiation 
and a scalar field evolve adiabatically without interaction among them. The evolution 
of a scalar field is governed by a quintessence, i.e. pq = WqPq, where Wq = — 1. 
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The quintessence field Q is a normalized scalar field \P. During a radiation era 
{W ~ 0) an equation of state for a scalar field ^ is different: pq, = ^pq, (similar 
to stiff matter equation of state). The scale factor R(t) evolves according to some 
elementary function built from hyperbolic function in radiation + quintessence scenario 
and according to complicated function formed from the first and the second elliptic 
integrals (after reversing of those functions). It expands and accelerates. The solution 
with radiation and quintessence has a restricted behaviour (Eq. (3.90-91)). The energy 
density of radiation is going to zero. Only a quintessence energy density is constant. 
The energy of a quintessence is a potential energy dominated. However, the solution 
cannot be continued up to the end for an interaction between a matter and a scalar 
field starts to be important. The effective gravitational constant remains constant 
during the period. Let us remind to the reader that during inflationary epochs the 
effective gravitational constant was constant. Only during a radiation era it was slightly 
changing. In order to keep the Universe to evolve forever it is necessary to find a 
solution with an interaction between a matter and a scalar field $ '. Using an appropriate 
ansatz for an evolution of an energy density of a matter we get such a solution. The 
scale factor R(t) is exponentially expanding. 

An energy density of a matter is going to zero and an energy density of a scalar 
field is approaching (exponentially) a constant. However, now a scalar field & forms 
a different form of a matter — a stiff matter with an equation of a state pq, ~ pq, 
(Wq, ~ 1). The effective gravitational constant is growing exponentially in time. In 
contradistinction with the previous period of an evolution the scalar field \P forms a 
matter with a kinetic energy dominance (i.e. -ftT-essence) . In the previous period we 
have to do with a potential energy dominance (i.e. quintessence). 

Let us give some remarks on a spatial curvature of a model of the Universe. In the 
first period (inflationary scenario) the spatial curvature has to be zero and any pertur- 
bations of initial conditions cannot change it. In next periods it is natural to suppose 
that it is not changed. Thus our model of a Universe is spatially flat. In the moment of 
change of a phase of an evolution we should apply a matching conditions for an energy 
density and for a scale factor (a radius of a model of a Universe). However, an evolu- 
tion in such moments undergoes phase transitions for a law of an evolution of a scalar 
field and a scale factor, which can be considered as a second order or even a first order 
phase transition. It is worth to notice that in formulae concerning Hubble constants, 
deceleration parameters, gravitational constants (effective) meet some parameters from 
unification theory of fundamental interactions and from cosmology. It seems that we 
are on a right track to give an account for large number hypothesis by P. Dirac. This 
demands of course some investigations especially in developing nonsymmetric Kaluza- 
Klein (Jordan-Thiry) theory. It is also interesting to mention that the scalar field ^ 
plays many roles in our theory. It works as an inflation field during an inflation epoch. 
(However Higgs' fields #~ play an important role in creating fluctuations spectrum.) It 
plays also a role of a quintessence and a K-essence, except its influence on an effective 
gravitational constant and on an effective scale of masses in unification of fundamental 
interactions. This field is massive getting masses from several mechanisms. During 
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inflation phases it acquires mass due to spontaneous symmetry breaking, different in 
both epochs. It gets a mass from a cosmological background. Thus except its cosmo- 
logical importance the fluctuations of a scalar field around its cosmological value could 
be detected (in principle) as massive scalar particles of a large mass. In this way an 
effective gravitational potential in nonrelativistic limit takes a form 

= G N M / _ « e -(^)\ eK(t - to) + Ar^ e - K (t-t ) (3 . 145) 
r V r J 6 

where k, A are given by Eqs (3.131-132), r is a range of massive scalar 5\P (a fluctuation 
of a scalar field If in a cosmological background). This could be detected as a tiny 
change of a Newton constant in gravitational law (e.g. the fifth force) on short distances. 
Let us notice that on short distances and on short time scales an effective gravitational 
potential reads: 

V{r) = ?^(l--e-W). (3.146) 
r V r / 

On large distances and short time scales 

_ + V 

r 6 

On short distances and large time scales 

y {r) = ^EMo ^ _ £ e -(^ gK(t _ H (3148) 

On intermediate distances and large time scales 

v{r) = GNMo eK{t - tol (3149) 
r 

The latest being a cause to collapse a cluster of galaxies into a black hole. It seems 
also possible to consider an effective coupling constant a e s s for a s enters and m e -?. 
However 



r is a radius of a manifold G/Gq (e.g. a radius of a sphere S ). If we rewrite Eq. (3.150) 

as 

-£ = -e"*, (3.151) 

it would be quite difficult to distinguish a drift of rrf^ from af drift. Only a quotient 
has a definite dependence on the scalar field if. Especially it is interesting to consider a 
drift of a em (a fine structure constant) reported from several sources. However, without 
extending our theory to include fermion (spinor) fields this is impossible. Thus it is 
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necessary to construct a nonsymmetric-geometric version with supersymmetry and 
supergravity including noncommutative (anticommutative) coordinates to settle the 
problem. The scalar field \P can form an infinite tower of massive scalar fields which 
could have important astrophysical and high-energy consequences. 

Let us consider Eqs (5.2.51) and (5.2.56) from the first point of Ref. [1], p. 318-319, 
and changing a base in the Lie algebra f) we get 

= St (3.152) 

and 

a\ = 51 (3.153) 
Simultaneously for G C H (q C f)) we get 

fk = Cl~ h (3.154a) 



ab ab 

= c ii ( 3 - i54b ) 



where a,b,c refer to the complement m in g (q = Qo+m), i, J, k to the subalgebra 0o 
of g. Let us suppose a symmetry requirement 



From Eq. (3.155) it follows that 

Eq. (2.5.51) looks like 
Eq. (2.5.55) 



[m,m]Cfl - (3.155) 

C L = fL = 0- (3-156) 
*£C& = #|C&. (3.157) 

H k = c «^H- c k-^ c i- ( 3 - 158 ) 

Let us come back to Eq. (3.1). 

For our inflation driving agent is a scalar field )P we do not carry any backreaction 
of Higgs' fields and we can consider some nonstandard scenarios of an evolution of 
those fields. Thus we can consider an evolution which is not a slow-roll evolution. Let 
us suppose that 

1 0*4)^0. (3.159) 

In this way we get from (3.159) 
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The next step in simplification is such that we collapse all the degrees of freedom of 
Higgs' fields into one 

(3.161) 



= SI (b 

m m ' 



In this way constraints (3.157) are satisfied trivially and Eq. (3.158) is as follows: 



m = c2r ($ 2 - 1 - $) 

ab ab V / 

In this way V = V (no constraints) and 



where 



and 



A 



V=^ (a> 2 - 1 - * s <pf 



/ I l )/i[" l ''ln[''']r ,c r<b fb „ah „bm t a \ 
lab[lg [ >9 [ l( ^rhh C d ~ b - C rhn9 9 



Kb = Kb ( 



1 



<2> 



1 r, r ,cm j d _i_ 7 n rhcjd o; „arh„mcsid 

l dc9rhb9 L cd + l cd9am9 ^~ b£ ~ ^cd9 9 <-^c' 

We can get the following identity 

ldc9 M 9 hp Li~ a C c M + l cd g^ b Lf £ C c M = 2l cd g^ h C%C c M . 
Using Eq. (3.160) we get the following equation for <P: 

d 2 <P e n *° 5V _ 

if the following constraints are satisfied 



J VW\d n x (g q jr + r*Wf ) = 2 ^<Pf- 



Thus from Eq. (3.168) one gets a first integral of motion 



<P 2 e rf V. B 



2 + 2r 



-V = — = const. 
2 



From Eq. (3.170) one gets 



I 



-A (a 2 H $ 2 _ as $ -if + B 

raj 



±{t-t ) 



(3.162) 

(3.163) 
(3.164) 

(3.165) 
(3.166) 

(3.167) 
(3.168) 
(3.169) 

(3.170) 



(3.171) 
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and finally 

/ , » =± J^L (t - to) (3.172) 

•Po 

where 



= 0!s Vi^ (3-173) 

<P = —<p. (3.174) 

a s 

In order to find an integral on the left hand side of Eq. (3.172) and its inverse 
function we use a uniformization theory of algebraic curves. Let 

f(<p) = (a-ip 2 + ^ + l)(a + ^ 2 -ip-l) = -ip i + 2i P 3 + ^ 2 -2^ + (a 2 - 1). (3.175) 

One easily notices that 

<p = (3.176) 

is a root of the polynomial f(<p) (a real root). 
In this way one gets 

, = <W-irw) +w (3 ' 177) 

where 

P(z) =P(z; 52 ,<7 3 ) (3.178) 
is a P Weierstrass function with invariants: 

92 = (3.179) 



3 

93 = <*- + 12 . (3.I8O) 
12 27 



One easily finds 



/'M = |(5 + 3a) v / 5 + 4o" (3.181) 
/"(^o) = -2(5 + 6a) (3.182) 



and 



z 



r[f(x)]-Ux. (3.183) 

J ion 
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In this way one gets 



*(*) = (5 + 3a)VCT IWOT 
W 8a s (P(z) + ^(5 + 6a)) 2a s V ; 



where 



z = ± ] J——(t-t ). (3.185) 

Let us come back to the potential V Eq. (3.163) and let us look for its critical 
points. One gets 

% = — (a^ 2 - a s $ - 1) {2a s <P - 1) (3.186) 



and from 



easily finds 



(3.187) 



<2>x = — (3.188a) 



<P 2 = — (3.188b) 

2a s 

<2>3 = (3.188c) 
2a s 



One obtains 



25 A 

m) = ^ (3.189a) 

V(<2> 2 ) = (3.189b) 
V(* 3 ) = 0. (3.189c) 

It is easy to see that #i is a maximum and <p2 and ^3 are minima of the potential V. 

However, in our simplified model of an evolution of Higgs' fields we have not a 
second (local) minimum. Thus in order to mimic a real situation we start an evolution 
of a Higgs field from the value 

<£ = 0. (3.190) 



In this case 



Thus 



* < : "4- « :U!H 



at 



^initial) = (3.192) 
^(tend) = ^^. (3.193) 
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Coming back to Eq. (3.172) we get 

1 + V5 



f ^ = \ ^—(tend -^initial) (3.194) 

J V(«-^ 2 + ^ + l)(a + ^-^-l) V r Vend mitiaU V ; 

Thus the amount of inflation obtained in our evolution is equal to 

1 + V5 
2 

dtp 



Supposing simply- 



one gets 



a > - (3.196) 



No = \h 



2" U V 2Aae n *° 

Vl , x \ — 1 V2 , , _ i \ (3.197) 



vO 



where 



or 



or 



cos2 "' = 5T^ (3 - 198) 
W » = f /O - (sTs) sm2 *r " 9 (3 - 200) 



¥>2 



In order to get 60-fold inflation 

N ~ 60 (3.202) 

we should play with parameters in our theory. Let us notice that in our simple model 
the constant a\ equals 

2 ( m A \ 2 f A\ 1 ( l p \\ 2 { A 



« - ™* ( = ^ [f U ) (3 ' 203) 
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(see Eqs (2.8-9a). 

Let us consider a slow-roll dynamic of Higgs' field in this simplified model. From 
Eq. (3.9) one gets 



d 
dt 



i}^ 9pn ~\~ I ^ 9hp) 



sv 



Using our simplifications from Eqs (3.152-162) and Eqs (3.143-149) one gets 



(3.204) 



d$ e"* 1 dV 



supposing a condition 



dt 12rH d$ 

&p 9pn 9np 



or 



d$ e n ^A 

-77 + 



a 



2 <P 2 - a a $-l) (2a s <2>- 1) = 0. 



dt 6a s rH 
Changing the dependent variable into 



ip = a s <P 



one gets 



where 



dip 
~dt 



+ b (if 2 - tp - 1) (2^9 - 1) = 



From Eq. (3.209) one finds 



e n* lA 

b = 1CJT > °- 
6rH 



and finally 



*(*) 



1 



a. 



11 15 — e~i 



2 2 y l - e -2 6 (*-*o) 
One easily notices that for (t — t ) —>■ 00 

1 / 1± y/E\ 



(3.205) 

(3.206) 
(3.207) 

(3.208) 
(3.209) 
(3.210) 

(3.211) 

(3.212) 



a. 



(3.213) 



It simply means that a slow-roll inflation never ends. 
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If we suppose that an inflation starts at (p = 0, i.e. for 

tinitiai=io-|rln2 (3.214) 

it is evident that to complete it we need an eternity. 

Summing up we get in our simplified model a finite amount of an inflation for a 
nonstandard dynamic of a Higgs field and an infinite amount of inflation for a slow-roll 
dynamic. The truth probably is in a middle. 

Probably it would be necessary to consider a full equation for a Higgs field in a 
simplified model 

^ - 3H ^- - 36# (a> 2 - a s <P - l) {2a s <P - 1) = (3.215) 

or 

^ - 3H ^ - 3ba s H {v 2 -V- 1) (2<p - 1) = 0. (3.215a) 
Changing an independent variable from t to r = H t one gets 

$-3^-3^(^-^-l)(2^-l) = 0. (3.215b) 



Let 



36a, 

= a (3.216) 



H 

and let us change a dependent variable <p into 

z = 2ip - 1 (3.217) 

V = ^ • (3-218) 

One gets 

ft - 3? - V - = 0. (3.219) 
Now let us take a special value for a = — ^ and let us change z into r: 

z = iy/5r. (3.220) 

One gets 

ft + 3^-2r 3 + 2r = 0. (3.221) 
This equation can be solved in general [10] 

r( T ) = iC ie T sn(Cie T + C 2 , -1). (3.222) 
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Thus 

z(t) = -V5de T sn(Cie r + C 2 , -1). (3.223) 

sn(u, —1) is a Jacobi elliptic function with a modulus K 2 = — 1, C\ and C 2 are 
constants. In this way 

${t) = ^—(l- V5C ie Hot sn(C ie Hnt + C 2 , -1)) ■ (3.224) 
2a s V / 

However, we are forced to put a = — y^. Let us remind that 

71^1 

'-sir 4 (3 ' 225) 

and we are supposing A > 0. However (see Eq. (3.164) for a definition of A), it is 
possible to consider A < only for a pleasure to play. Thus we use formula (3.224) to 
consider a dynamics of a Higgs field. Let us start an inflation for t = with <P = 0. 



One finds 



(i.e. initial = 0). 
Let 



<2>(0) = 0. (3.226) 



1 sn(Ci + C 2 ,-l) (3.227) 



n, y 2 

(a true minimum — a "true" vacuum). One gets 

1 



*(*end) = — ] (3.228) 



= -Cie Hntend sn (Cie Hntend + C 2 , -l) . (3.229) 



2 

Let us calculate the derivative of @(t). One gets 



f =-^C.^(-(Ci«*« + ft,-l) 

+ cn (de* 1 + C 2 , -1) • dn (Cie ff » f + C 2 , -l)) . 



(3.230) 



Let us suppose a slow movement of such that 



— (0)=0. (3.231) 
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From (3.231) one gets 

sc(Ci + C 2 , -1) = - dn(d + C 2 , -1). (3.232) 

This gives us an equation for a sum of integration constants. Thus we can get from 

Eq. (3.232) and Eq. (3.227) integration constants C\ and C 2 and from Eq. (3.229) t en d, 
which can give us an amount of inflation 

N = H t end . (3.233) 

In some sense Eq. (3.229) is an equation for an amount of inflation 

I = -C ie N ° sn (c ie Ra + C 2 , -l) . (3.234) 
Using Eq. (3.232) and Eq. (3.227) one easily gets 



d = -y 5 ^ ~ -0.5256. (3.235) 

One gets 

d + C 2 = -= (F(90°\45°) - F(32°\45 )) = 0.9431 (3.236) 
v2 

where F is an elliptic integral of the first kind and 



arcsin j \j ^ ) ~ arcsin (0.52.")(i ) ~ . (:i.237) 



Let us come back to Eq. (3.234) and let us denote 

e Ro = x . (3.238) 

One gets using (3.236) and (3.235) 

C 2 ~ 1.4687 (3.239) 

and finally 

1.4866xo - 5.0354 = F (^/45°) , (3.240) 

where 

1.9026 . A . 

cos 09 = . 3.241 

x 
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From Eq. (3.240) one finds 



2 ' 8284 5.0354 / ■ <"> (3.242) 



We come back to this equation later. 

Thus we get the following results. The Higgs field evolves from a "false" vacuum 
value (<P = 0) to the "true" vacuum value (<P = <Z>o) completing a second order phase 
transition. The potential of selfinteracting \P field changes and a new equilibrium <Po is 
attained. However, the field \P evolves from $b to new value a little different from $0 
and afterwards a radiation era starts. The change of the value of a scalar field \P is 
small in terms of a variable y (see Eq. (2.69) for a definition) only from to 1 (n is 
a natural number n = dim//, and H is at least G2, dimG2 = 14). Let us consider an 
evolution of a field ^ in this epoch. From Eq. (2.59) we get 



ip + u%<p = (3.243) 



where 



V = & + <p (3.245) 

and we linearize Eq. (2.59) around <Pb neglecting a term with a first derivative of ip. 
Thus a scalar field \P undergoes small oscillations around the equilibrium These 
oscillations cannot be too long for a field should change from 



* = ta (VCTj (3 ' 246) 

to 

W = ln(Jtpj. (3.247) 

In terms of a field tp from (p = to <p = ^ ln(l + — ) ~ ^. Thus we have 

<p = <Po sm(u t). (3.248) 
For t = 0, (f = the time Z\t to go from ip = to 92 = ^ is simply equal to 

arcsin(^— ) 



At = 



u 2u 
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< . (3.249) 



Thus the full amount of an inflation in this epoch is equal to 



N x = AtHi = — 



arcsm 



ip n 



< 



vrm pl 



2 ^o 2V6M(n + 2) 



(3.250) 



In our notation 



t 1 =t r + At 



(3.251) 



(see Eq. (2.89)). 

After a time t\ the Universe goes to the phase of radiation domination with a strong 
interaction between a radiation and a scalar field ^ up to a minimal temperature where 
an ordinary (a dust) matter appears. This matter evolves afterwards independently of 
a radiation and of a scalar field $ '. The scalar field now is evolving as a quintessence. 

In order to get some comparison let us consider an evolution of a scalar field & 
during the second de Sitter phase in a slow-roll approximation. Using Eq. (2.59) one 
gets 



dy 
dt 



(3.252) 



or 



where 



\A - V 2 dy 



5 = 2(n + 2)| 7 | 



n 



n + 2 



-B{t-t 



(3.253) 



(3.254) 



and 




(3.255) 



Let 



J y ( ^V-(^)) J 



for 



\A - y 2 dy 

(y 2 - (t^)) 



(3.256) 



n 



<y<i 



n + 2 
n > dimG2 



14. 



(3.257) 
(3.258) 
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For the integral / one finds: 
I = — arcsin(y) + 




(3.259) 



In order to find an amount of an inflation let us find limits of I for y = 1 and y = ^J^^. 
One finds 



lim / 



IT 

-2 + 



2 

n 



2 

n 



(n + 2 
In [ 2 (l - 



n+2 J J 



(n + 2) 



lim / 

, / n 
V " + 2 



-oo. 



(3.260) 
(3.261) 



Thus we see that a slow-roll approximation offers an infinite in time evolution of a 



field \P from y ^2 to ^ ^ ^ s s i mnar to an evolution of a Higgs field in some slow-roll 
approximation schemes. The inflation never ends. Let us come back to the Eq. (2.110) 
and let us calculate a Hubble constant and a deceleration parameter for this model 



_ R _ ( 3M(3 n (n- 1) 
~ R ~ V 16|7|"+ 4 



1 + 



16|7|"+ 4 
3M/?"(n- 



T)(*-^) 2 ) 



g = -- 



3Mf(n-l) 



(t-ti)" 2 . 



(3.262) 
(3.263) 



R 2 16| 7 |"+ 4 

Let us compare Eq. (3.262) with the similar equation for radiation dominated 
Universe in General Relativity 



(t-h 



and let us calculate a speed up factor 

H 3M(3 n (n-l) 



^1 = 



(t-h 



H GR 16| 7 |"+ 4 A , 16|7l"+ 4 (t t )2 \ 



(3.264) 



(3.265) 
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and let us do the same for a model Eq. (3.58) (i.e. radiation dominated with a quint- 
essence). Using Eq. (3.59) one gets 



A 2 = — = V^ m P iM*-*i) ( 3 266) 

#GR B tgh ( 2 -^m vX (t - ' 



Let us notice that for a (t — ti) ~ 



A 1 ~ (3.267) 



and 



9(n - l) 2 M 2 /3 2n 



lim = or^, i, ,« • (3.268) 
t^oo 256|7| 2n + 8 



Thus at early stages model Eq. (2.110) is slower than that in General Relativity. 
For large time the behaviour depends on details of the theory. In the case of model 
Eq. (3.58) we have 



and 



lim A 2 = \ (3.269) 

t — >t\ 



lim A 2 = +oo. (3.270) 

t— >oo 



Thus in the first case A is monotonically going from to 9 ^256|7|^"+^ — an< ^ m ^ ne 
second case from ^ to infinity. According to modern ideas an expansion rate in early 
Universe has an important influence on a production of light elements, the so called 
primordial abundance of light elements (see Ref. [11]). If at a beginning of primordial 
nucleosynthesis the Universe expansion rate is slower than in GR, then we have 4 He 
underproduction. This can be balanced by considering a larger ratio of a number 
density of barions to number density of photons (remember the barion number is a 
conserved quantity) 

7] = — (3.271) 

where ub, n 7 — barion and photon density numbers for a high redshift z = 10 10 . In con- 
trast to the latest if an expansion rate became faster than in GR during nucleosynthesis 
process, those bigger rj (traditionally one uses the so called) 7710 

7710 = 10 10 ?? (3.272) 

do not result in excessive burning of deuterium because this happens in a shorter time. 
The standard model of big-bang nucleosynthesis (SBBN) demands 

3 < mo < 5.6. (3.273) 
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It seems in the light of observational data from cosmic microwave background (CMB) 
(BOOMERANG and MAXIMA) and from the Lyman a-forest that r] W = 8.8 ± 1.4, 
significantly higher than the SBBM value (for CMB) and 8.2 ±2 or 12.3 ±2 for Lyman 
a-forest. Thus the second model could help in principle to explain the data without 
modification of nuclear reaction rates due to neutrino degeneracy or introducing new 
decaying particles. As the Universe expands, cools and becomes more dilute, the 
nuclear reactions cease to create and destroy nuclei. The abundances of the light 
nuclei formed during this epoch are determined by the competition between a time 
available (an expansion rate) and a density of reactans: neutrons and protons. 

The abundances of D, 3 He, 7 Li are limited by an expansion rate and are deter- 
mined by the competition between the nuclear production/destruction rates and an 
(universal) expansion of the Universe. The SBBN theory is based on a flat radiation 
dominated Universe model in General Relativity and found in a laboratory nuclear 
reaction rates. Thus it is strongly constrained. Any significant discrepancy between 
observed and calculated value of 7710 (known as baryometry) could be very dangerous. 
Thus changing the ratio of an universal expansion relative to the model of General 
Relativity can give some margin in a primordial alchemy. 

In our theory after two inflationary phases we have in principle two radiation 
dominated phases described by Eq. (2.110) and Eq. (3.58) with speed up factors (3.265) 
and (3.266). However, the important point is to match those models. We get 



R(t d ) 



R 1 exp {H t r + H^h - t r )) ( 1 + 77. 



2 

min 



- I ah I - m pl hi - t + — w r, min 



(3.274) 



and 



B = Pr{td)R (t d ) = pr *1 + ; |2± 4 ??min \R *1 + , , n±4 

V 2 | 7 | 2 J \ 2| 7 | 2 



Pmin Rl exp (3H t r + SHtih - t r )) 1 



(n + 1) 



^?min 



where p m i n is given by 



rn 



ont 

2 P n , „ \2(n+l) 



Pl| 



(l+Tfeninr"" 1 " 'WkCl + Tfcnin), 



(3.275) 



(3.276) 



^4(7/) is given by formula (2.90) and ?7 min by formula (2.97). 

The time t r = t enc j in any inflation model for an evolution of the Higgs field in the 
first de Sitter phase. The time t\ 



t r + At can be calculated from Eq. (3.249) 



arcsin(^— ' 
At = 



(3.277) 
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and is bounded by 



*<- = -*-M l^rA— Y(A) ■ (3-278) 



From Eq. (3.274) we can get the constant t 

~M, 
2M 



*0 = *1 + ~ .+4 ?7min - T (3.279) 



2 



where 



i?arsh 
T = 



flf^exp (2fl"o«r + 2ffl(tl - * r )) (l + ^Vlin) 



2\/3 



(3.280) 



and B is given by Eq. (3.275). 

Let us notice that for t = td the Universe undergoes a phase transition due to a 
change in an equation of state for a scalar field \P from pq, = 3p# in a first radiation 
dominated phase to p& = —py in a second radiation dominated phase (a quintessence 
phase). 

A matter which appears in the second phase does not play any role in an evolution 
of the Universe. For t = td we have a discontinuity in Hubble constants (parameters) 
for both phases 

H{t d ) = • Vm _ in - (3.281) 

X 1 + (n-l^min J 

in the first phase and 

H{td) = V3m pl (^jp) ctgh (^^m^ (3.282) 

for the second one, 

H(t d )^H(t d ) (3.283) 

1 A (M,\- l^l t+lnt 
-Aco^oj — sn. 

2 p 2 (n + 2) 2 

Let us notice that a speed up factor A\ changes from 



PQ = tjAcoOPo) = H¥7I I ^¥TT - ( 3 - 284 ) 



^i(ti) = (3.285) 

to 

^i(*d) = ^ • 7 r- (3.286) 



{} + (n-l^min) 
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and an expansion rate is slower than in General Relativity. 

The speed up factor A2 has a more complicated behaviour than A\. It is natural 
to expect a rapid speed up of a rate of expansion resulting in non-equilibrium nuclear 
synthesis of light elements. Thus the presented scenario offers interesting possibilities 
for an evolution of early Universe. Recently some papers have appeared on scalar- 
tensor theories of gravitation exploiting the idea of a speed up factor in primordial 
nucleosynthesis [12]. 

Finally let us come back to the model (3.64), (3.67). We cannot invert analytically 
the formula (3.67). Moreover taking under consideration Eqs (3.85-88) and making 
some natural simplifications we come to the following approximative formulae for R(t). 

For x < 1.1969 (Eq. (3.86)) 

/ A 

R(t) = J — (0.44 + 0.0857V) (3.287) 
V PQ 

where 



N = exp ( 0.374^^(t - t ) + 157.93 ) . (3.288) 

V m Pi ) 

One can calculate a Hubble parameter and a deceleration parameter and gets: 

H — — — M 3A7N (3.289) 
R mpi 44 + 8.5iV v ; 

RR 5.21 

-q = = -^7- + 1- 3.290) 

We have H > 0, and q < 0. Thus the model expands and accelerates. First of all we 
consider 

= A e «Q 
p m R3 e 

and we take boundary for R. In this way one gets 

p m = a l -p Q e aQ , i = l,2,3,4, (3.291) 

ai = 2.015 
a 2 = 0.034 

3.292 

a 3 = 1060 V ; 

04 = 8 

If we reconsider a contemporary gravitational constant Gn as GAre~( n+2 ^ n we would 
get interesting ratios. Thus we get 

— = 0.034 -7- 2.015 or 8 H- 1060 (3.293) 
PQ 
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which for the first is an excellent agreement with recent data concerning an acceleration 
of the Universe. Solving Eq. (3.288) for TV if R = 1-33^J gives us 

N = 10.47 (3.294) 



H = • 0.25 (3.295) 
m p] 

-q = 1.5. (3.296) 

Using Eq. (3.288) and Eq. (3.295) we can estimate a time of our contemporary 
epoch 

t = -i-10 15 yr-to (3.297) 

where h is a dimensionless Hubble parameter, H = h- H (we take for a contemporary 
Hubble parameter Hq = 100^^), 0.7 < h < 1, which seems to be too much. 
We can also estimate a density of a quintessence 

H 2 ks £ 

po = — ■ 0.571 = h 2 ■ 0.08 • 10 _23 -| = 8h 2 ■ 10 -29 -^- . (3.298) 
Gn nr 3 cnr 3 

Let us come back to the Eqs (3.177) and (3.184). In Eq. (3.184) a Higgs field is 
given by a P Weierstrass function. This function can be expressed by Jacobi elliptic 
function 

P(z) = e 3 + (ei - e 2 ) ns 2 [z{e 1 - e 2 )*) (3.299) 



with 



K 2 = (3.300) 

ei - e 3 



where ei, e 2 , e 3 are roots of the polynomial 

Ax 3 -g 2 x-g 3 . (3.301) 
Thus we should solve a cubic equation 

1 „ on 1 /a 2 10 N 



3 



i2 ( 5 - 3 ° 2 ) - 12 T + y) =a (3 - 302) 



We want Eq. (3.302) to have all real roots. Thus we need 



12 

and 



p = --3- (5-3a 2 ) < (3.303) 



„3 „2 

D = tj + T < (3 - 304) 
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where 

1 fa 2 10 

Both conditions (3.303) and (3.304) are satisfied if 

< a < 0.3115. 



In this case we get 



where 



From (3.307-309) one gets 



ei = —y/5 — 3a 2 cos ^ 
o o 

e2 = - V 5 — 6CL Z cos — - — 

= - v 5 — 3a^ cos 

3 3 



(9a 2 + 40) 

COS(f - 



12(5 - 3a 2 ) 3 / 2 



sin ^ 
if 2 = 3 



sin ^ — 



Thus 



where 



ei > e2 > e3 



m ^ (5 + 3a)^/5T4a" (1 + V5T4a) 



JAe n *° (5Ae n *° - 3a 2 r) / <p + vr . 



/5 9 . V 3 + 7r 
-Ai = ei — 62 = y 77 — a z sin — - — 



and 



1 . , / 5 „ <f + 47T 1 . 

^2 = e 3 + — (5 + 6a) = J - - a 2 cos + —(5 + 



COS </3 = 



12 (5Ae n *° - 3a 2 s Brf /2 



a 2 Br 
< -rS-z- < 0.3115 
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and a is given by Eq. (3.173). 

The function ns(u, K) is an elliptic Jacobi function with a modulus K given by 
Eq. (3.311) 

n S (u,K)= 1 (3.319) 
sn(u, K ) 

sn(u, K) = sinam(u, K). (3.320) 

In order to justify our treatment of Eqs (3.130) and (3.133) we consider the full field 
equations 

R^u ~ -^9tiuR = 2Au^u„ - Ag^ u (3.321) 

(where A is given by Eq. (3.136)). 

However, in this case we consider A arbitrary. Let us consider static and spherically 
symmetric metric 

ds 2 = e v dt 2 - e A dr 2 - r 2 (d6 2 + sin 2 6 dip 2 ) (3.322) 

where v = v{r) and A = A(r). 

From Eqs (3.321-322) one gets 

r 2 A = e - x {rv' + 1) - 1 (3.323) 

r 2 A = e" A (rA' - 1) + 1 (3.324) 

A.' = -Av' (3.325) 



where ' means derivation with respect to r. Summing up (3.323) and (3.324) one gets 

e" A (A' +v') = 2Ar. (3.326) 
Moreover in our case A is very small. Thus we get approximately 

^ ~ (3.327) 
dr 

and 

A' -v'. (3.328) 

In this way we go to the solution (3.137) for a small A which is our case. In this way 
A = — v and 

B(r) = e v = e~ x . (3.329) 

Let us check a consistency of our solution. In order to do this we consider Eqs 
(3.323-325) in full. One gets from Eq. (3.325) 

A = e-^e~ v . (3.330) 
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For A is very small we should suppose that fio is positive and very large (fj,o — > oo). 
From Eqs (3.323-324) one gets 

4-(* + v) =2e-^ n ■r-e (A -" ) . (3.331) 
ar 

Supposing that 

A(r ) + w(r ) = (3.332) 

for an established r > (r is greater than a Schwarzschild radius of the mass M 
from the solution (3.137)). Taking sufficiently big no we get approximately 

^(A + w)~ (3.333) 

and of course 

A -v. (3.334) 
Let us come back to the Eq. (3.242) and consider it for 

(p = kir + ^-e, fc = 0,±l,±2, ... , (3.335) 

where e is considered to be small. One gets 



2 - 82M ^ - 5.0354 



sine 



F 

2 fc 

/ ^ + 2fcifQ), (3.336) 



where 



d() 

K{\)= j = 1.8541 (3.337) 

'l - i sin 2 



V 2 



is a full elliptic integral of the first kind for a modulus equal ^ . For a small e we can 
write 

sine~e (3.338) 

/ = = if (|) — v^e. (3.339) 

Thus one gets for k = 21 

2.8284 r- 

5.0354 = -V2e + (4/ + 1)1.8541. (3.340) 

£ 
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It is easy to notice that we can realize a 60-fold inflation in our model for 

1.9026 , , , 

x = ; ; 3.341 
|e| 

can be arbitrary big for sufficiently big I. 

Let us take I = 25 ■ 10 24 . In this case we have for e an equation 

e 2 - ((41 + 1)1.3110 + 5.0354)e + 2 = (3.342) 

e 2 - (10 26 + 6.0354) e + 2 = 0, (3.342a) 
e 2 - 10 26 e + 2 = 0, (3.342b) 



or 



i.e. 



and finally 

|e| ~ 10" 26 (3.343) 

x ^ 1.9026 • 10 26 (3.344) 

lnx ~ 60.51 (3.345) 

which gives us a 60-fold inflation. 

Moreover the Eq. (3.242) has an infinite number of solutions for < tp < ^, 

ip = kir + ip, 

- 5.0354 = I , d ° + 7.5164 • / (3.346) 
c ^ { ^l-^in 2 e 

with 

1 9026 

x = - , 1 = 0,1,2,3,... , k = 21. (3.347) 

cos tp 

One can find roots of Eq. (3.346) for 
/ = 5 

x = 16.83 
lnx = 2.82, 

/ = 50 

x = 0.24 • 10 3 
lnx = 5.48 

and lnxo for / = 5 • 10 24 . In the last case one finds using 70-digit arithmetics 



(3.348a) 
(3.348b) 



In x ^ 58.479... (3.349) 
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which gives us an almost 60-fold inflation. In general an amount of inflation 

iV =lna;o (3.350) 

is a function of / = 0, 1, 2, . . . and probably can be connected to the Dirac large number 
hypothesis. 

Finally let us take I = 10 n where n > 10. In this case one can solve Eq. (3.342a) 
and get 

2 

|p| ~ C3 351) 

1 1 [(4 -10™ + 1)1.3110 + 5.0354] ' v ■ J 

Thus for x we find 

x ~ 5 • 10" (3.352) 

and 

N (n) ~ 1.60 + 2.30n. (3.353) 
In this way, for large I, N is a linear function of a logarithm of I, 



iVo(10) ~ 24.6, iV (20) ~ 47.6, 7V (24) ~ 56.8, 
iV (25) ~ 59.1, iV (26) ~ 61.4 



(3.354) 



or 

7V (log 10 Z) = 1.6 + In/. (3.355) 

It is easy to see that Eq. (3.355) is an excellent approximation even for I = 50. 

Let us come back to the Eq. (3.29) in order to find a power spectrum for our 
simplified model of inflation. Using Eqs (3.161), (3.230), (3.235) and (3.239) one gets 
after some algebra 

ft<*)* 2.8944 (g) 2 „ ia ;/(^) (3.356) 

where 

f{x) = x~ 2 (sn(1.4687 - 0.5256x, -1) 

+ cn(1.4687 - 0.5256x, -1) dn(1.4687 - 0.5256x, -1)J . 
Using some relations among elliptic functions one gets 
, x 1 2dn 4 (u, i) 

f{x) = 2 ( 3 - 358 ) 

x2 (sn( U ,i)dn(n,i) + v / 2cn( U ,i)) 

where 

u = 2.0770 - 0.7433a;. (3.359) 



(3.357) 
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Let us consider a more general situation for the Eq. (3.231), i.e. 



u (0) = h (3.360) 



where h ^ 0. In this way one gets 



h = ~y~ C 1 F (sn(C 1 + C 2 , -1) + cn(d + C 2 , -1) dn(d + C 2 , -1)). (3.361) 

a s 

Using Eq. (3.227) one gets 



Co 



V2 



~ 75 Jo 



\/\/(l + ^) 4 + 4-(l + ^) 

cos y — g 



m 2 

dip 



(3.362) 



1 — i sin 2 99 



and 



'10 

|d I = = . (3.363) 



yV( 1+ *) 4+4 -( 1+ *) 2 



Let us come back to the Eq. (3.229) to find an amount of inflation for C 2 and C\ 
given by Eqs (3.362-363). One gets 

arccos (^7^r) 

^2C 2 -K{\)= [ . d ° -V2C ie N \ (3.364) 

or using a natural substitution 



'1-isin^ 



1 

COSV9= 2C^' (3 ' 365) 



J Ji-i s ^e 2cos( ^ 



(3.366) 

V 2 

Taking as usual 

^ = 2/vr + f-e, 1 = 0,1,2,..., (3.367) 
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one gets for small e 



~ V2{Al + 2)K{\) - 2C 2 (4/ + 2)V2K(\) 
(for large I) and eventually 

Let us calculate ^ for i = t* = ln( R ^ o ). One gets 



■ft* 

tit v y 2a. 



+m ( c '(i^) +c -- 1 ) dn ( c -(i^) +c -- 1 ))' 

Thus we can write down a Pr(K) function. 



2vr y 5C2 ' V^o^o 



where 



f(x) = ^ (sn(Cis + C 2 , -1) + cn(dx + C 2 , -1) dn(dx + C 2 , -1)) 
Using some relation among elliptic functions one finds 

where 

1 



(sd(«, \) + >/2 cd(u, |) nd(n, \)) 
and 

it = V2~dx + \/2C2 . 
Moreover we can reparametrize (3.374-375) in the following way 

arccos ( ^ ) 

= V2C 1 x + K{\)-C 1 V2- I . dlf 



u 



o V 1 - f sm2 V 
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where 



For large C\ one gets 



\C X \ > VE. (3.378) 



u^V2d{x-l) (3.379) 



a s h C\ 



(3.380) 



Using (3.369), 



^ <*-i>. (3.3 81 ; 



If we take large C\ (it means, a large h) 

h = %^ (3.382) 

and simultaneously sufficiently large I we can achieve a 60-fold inflation with a function 
Pr(K) given by the formula (3.374). Large C\ means here C\ ~ 100, large I means 
I ~ 10 25 . 

Let us calculate the spectral index for our Pr{K) function, i.e. 

^ } ~ 1= dlnitT ■ (3 - 383) 

fin 1 9 9^ vgcd^j) - sd 3 (n, \) 

n s (K)-\ = -2-2C 1 x — (3.384) 

sd(tt, ±) + v2 cd(u, i) nd(u, ±) 

where 

« = V2(C\x + C 2 ) = -5^- (C\K + C 2 i?o#o) 

itO-Ho 

X = 



One gets 



(3.385) 



dn s (K) 

A very interesting characteristic of Pr{K) is also s — . One gets 
dn s (K) v/2cd(u, i) - sd 3 (u, ±) 



dlnK sd(«,i) + v^cd(u,i)nd(u,i 

-2C 2 x 2 sd 2 (n, i) -2Cfz 



2 ..2 .,2, , ,2 .2 (v / 2cd( U ,i)-sd 3 (n, i)) 



2 



(3.386) 



(sd(u, |) + >/2cd(u, |) nd(u, 
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where u and x are given by Eq. (3.385). 

Let us take for a trial C\ = —0.5256 and C2 = 1.4687. In this case one finds 

n s (if)-l = -2 + 1.0512x — ^^ / ' 2 ' (3.387) 

sd(«, 5) + a/2 cd(u, 3) nd(u, |) 

dn »(*0 _ x 0512a; V2cd(u, \) - sd 3 (u, 1) 



din if " sd(u,i) + ^/2cd(u,i)nd(n,i) 

" 1 <« 1 x2 (3-388) 

9 o / 1 \ o (\/2cd(u i) -Sd 3 (u I)) 

- 0.5525x 2 sd 2 (u, ±) - 0.5525x 2 ^ 2 

(sd(«, |) + >/2 cd(u, \) nd(n, f )) 

u = 2.0771 - 0.7433x 

K (3.388a) 

RqHq 

The interesting point is to find n s (K) ~ 1 (a flat power spectrum). One gets 

sd(u, ±) + V2 cd(u, |) nd(u, §) = Cix (sd 3 (u, §) - v^cd(u, (3.389) 

and 

sd(u, i) + y/2 cd(u, |) nd(u, ±) / 0. (3.390) 
Using (3.389-390) one gets 

^ = -2C 1 Vsd 2 ( U ,i) (3.391) 

if Eqs (3.389-390) are satisfied. 

In the case of special C\ and C2 one gets 

-0.5525x 2 sd 2 (u ±) (3.392) 



dlnK ' v ' 2/ 

where x, u are given by Eq. (3.388a). 

Let us reparametrize the Eq. (3.389). One gets 

sd(u, \) + V2cd(u, \) nd(u, \) = U ~^ 2 (sd 3 (u, \) - >/2cd(u, §)) . (3.393) 

For sufficiently big u one gets (the equation has infinite number of roots) 

u = u 1 + 2lK(\), I = 0,±1,±2, ... (3.394) 
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where u\ satisfies the equation 



i 

sn(«i,i) = ^ ^1 + ^5^65 + ^20357- \Av/20357 - 65^ " ~ 0.65... (3.395) 



Moreover for x = V K „ > we have the condition 



ui + 2/K(|) - y^Ca 



One finds 



> 0. (3.396) 



sd 2 (m,i) = 0.53... (3.397) 

Thus 

_(„ 1 + 2/^))=. -0.53^ ^ j , (3 ^ 98) 

I = 0,±1,±2,... 

One gets 

arcsin0.65 ~ 40°. 54 (3.399) 

40°. 54 

/df) 
-_ = F (40° .54/45°) ^ 0.73 (3.400) 
^1-isin 2 ^ 

Taking special value of C\ and C 2 one gets 

x 1.81 -4.68/ (3.401) 

where I = 0, ±1,±2, ... . For x > it is easy to see that / should be nonpositive. 
Practically / should be a negative integer 

I < -3 (3.402) 

and 

-^-^ =■ -0.29(1.81 - 4.68/) 2 . (3.403) 

Thus for 

Ki = R H (1.81 - 4.68/) 

we get 

= 1. (3.404) 
The important range of In K , A In if is about 10. 
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Thus 

n s (K) ~ 1 ± 2.9(1.81 - 4.68/) 2 . (3.405) 

Taking I = — 3 one gets 

-300 < n s (K) < 300 (3.406) 

and 

P R (K) ~ if 1_7l "(- K ") 

is not flat in the range considered. 

Moreover we can improve the results considering Eq. (3.398) for large C\. For 
large C\ one gets 

£r = -1. (3-407) 

Thus we find 

_^£_ ~ _i. 06 ( Ul + 2Zif (I) + d) 2 . (3.408) 
Taking large value of Ci in such a way that 

Ci = -21K{\) - Ul +e (3.409) 

where e is a small number, 

e ~ HT", 

one gets 

-1.06 -lO" 2 ". (3.410) 



din if 

The last condition means that we should take 



h ^ flojC^J = _Ho_ / .73-2/if(|)-10-") (3.411) 

V5a s V5r' 



where Z is an integer, / ~ —100. 

Thus for some special values of h we can get arbitrarily small — — a — which means 

a m A 

we can achieve a flat power spectrum in the range considered (Z\lnif ~ 10), 



n 



,{K) = l±10- (2n_1) , (3.412) 



n arbitrarily big. 

Let us consider the value of K coresponding to our value of n s (K) = 1. One gets 

0.73 + 21K(\) - V2C 2 

x = . 3.413 

V2C, 
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Using our assumption on a large C\ and Eq. (3.409) one finds 



1 



V2J 21KU) 



(3.414) 



or (for e is small and I quite big) 

x ~ 0.293 (3.415) 

and 

K ~ 0.3R H Q (3.416) 
with the range Zilni^T ~ 10. It means that 

0.3e" 10 < — ^— < 0.3 • e 10 (3.417) 

which gives us a full cosmologically interesting region 

K 

10" 4 < — — < 10 4 . (3.418) 
RqH 



Conclusions 

In the paper we consider some cosmological consequences of the Nonsymmetric 
Kaluza-Klein (Jordan-Thiry) Theory. Especially we use the scalar field \& appearing 
there in order to get cosmological models with a quintessence and phase transitions. 
We consider a dynamics of Higgs' fields with various approximations and models of 
inflation. 

Eventually we develop a toy model of this dynamics to obtain an amount of infla- 
tion and Pr{K) function (a spectral function for fluctuations). We calculate a spectral 

index n s (K) and — — s — for this model, 
am A 
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